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Abstract 

In the present model, a unified picture of cosmology from early in- 
fiation to late acceleration is obtained from f{R)— gravity with non- 
linear terms and of scalar curvature R. It is discussed that 
elementary particles and radiation are produced during early infia- 
tion. The emitted radiation thermalizes the universe to a very high 
temperature ~ lO^^GeV. The exit of the universe from infiationary 
phase is followed by deceleration due to radiation-dominance heralding 
the standard cosmology with background radiation having the initial 
temperature ~ lO^^GeV. It is found that dark matter is induced by 
curvature and baryonic matter is produced during infiation. Radiation- 
dominated phase is followed by deceleration due to matter-dominance. 
The curvature-induced phantom dark energy dominates at the red-shift 
z = 0.303 causing the late acceleration. It is found that the universe 
will collapse in future. Further, it is investigated that back-reaction of 

quantum particles produced near the collapse time can avoid cosmic 
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collapse and the universe will escape to revival of the state of early 
universe. 

Key words : f{R)— cosmology, early inflation, late acceleration, pro- 
duction of elementary particles, cosmic collapse and its quantum avoid- 
ance. 



1. Introduction 

Results of astrophysical observations, during last ten years, have con- 
clusive evidence for the late cosmic acceleration [1, 2]. It is caused by 
some exotic fluid with negative pressure p < — l/3pde, where pde > is 
the density for the dark energy fluid. This condition violates the cosmic 
strong energy condition (SEC) or weak energy condition (WEC). 

Many field-theoretical and hydro-dynamical models of DE were pro- 
posed in the recent past to explain the challenging observation of late 
acceleration. A detailed review of these models is available in [3]. In this 
race, curvature was also used as an important candidate and various 
models were proposed, where non-linear curvature terms were consid- 
ered as gravitational alternatives of DE. These f{R)— DE models are 
reviewed in [5]. Recently, these models specially the model with non- 
hnear terms i?" and R"'^ [6] received criticism in [7] on the ground that 
these do not produce matter in the late universe needed for formation 
of large scale structure in the universe. In another review, Nojiri and 
Odintsov have discussed dark matter in the late universe [8] responding 
to this criticism [7]. 

Here, we have a deviation in approach to obtain DE from curvature. 
In this approach, which was used earlier in refs.[9, 10, 11, 12, 13, 14]. 
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It is contrary to the approach in refs.[5, 7, 8]. Here, non-Unear curva- 
ture terms are not assumed as DE Hke refs.[5, 7, 8]. Rather, DE terms 
emerge from curvature spontaneously in this approach. In refs.[5, 7, 8], 
gravitational equations are derived from the action having Einstein - 
Hilbert term and non-linear curvature terms. Terms in gravitational 
equations, due to non -linear curvature terms in the action, arc recog- 
nized as DE terms. In the present paper as well as refs.[9, 10, 11, 12, 
13, 14], trace of /(i?)— gravitational equations are obtained yielding an 
equation for scalar curvature R. In the homogeneous space-time, the 
equation for R reduces to the second-order equation for the scale factor. 
First integral of this differential equation yields the Friedmann equa- 
tion giving dynamics of the universe. Here, it is found that DE terms 
emerge in the Friedmann equation so derived. Thus, in the present 
model, DE terms are induced by linear as well as non-linear terms of 
scalar curvature, whereas only non-linear terms of curvature contribute 
to the dark energy in [5, 7, 8]. 

Although the present model is based on f{R)— gravity, it is differ- 
ent from models in [5, 7, 8] in two ways. The main difference is the 
difference in approach mentioned above. The other difference is the 
advantage of getting matter term from the gravitational sector. More- 
over, there are many other interesting features of the present model. In 
the present set-up, non-linear terms of curvature are taken as R^ and 
7^(2+'-) along with Einstein-Hilbert term. This paper is an extension 
of the work [13]. Here many other aspects of cosmology such as cre- 
ation of SM particles in the early universe, collapse of the universe in 
future and possible avoidance of future cosmic collapse using quantum 
gravity are addressed to. In what follows, it is found that quintessence 
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dark energy and dark radiation are induced by curvature in the early 
universe as well as dark matter and phantom dark energy are induced 
in the late universe if r = 3. Here, the topology of the cosmological 
model is given by homogeneous and flat FRW universe. 

The present model yields an interesting cosmological picture from 
the early universe to the future universe, which is given as follows. 
Here, investigations begin from the Planck scale, which is the funda- 
mental scale. It is found that the early universe inflates for a short 
period, driven by curvature-induced quintessence scalar behaving as 
an infiaton. During this period, elementary particles are produced and 
lot of energy is released . The emitted energy thermalizes the universe 
rapidly upto the temperature ~ lO^^GeV. As a consequence, produced 
particles due to the decay of curvature-induced quintessence attain the 
thermal equilibrium with radiation. There are two sources of radia- 
tion (i) emitted radiation during inflation and (ii) curvature-induced 
radiation recognized as dark radiation. The density of the emitted ra- 
diation during inflation dominates, so the universe is driven by emitted 
radiation and decelerates after exit from inflation. Thus, in the present 
scenario, the standard cosmology is recovered at this epoch . The emit- 
ted radiation, during inflation, is identified with the cosmic background 
radiation with very high initial temperature ~ lO^^GeV. 

The produced elementary particles, during inflation, undergo various 
processes of the standard cosmology such as nucleosynthesis, baryosynthesis 
and hydrogen-recombination, which are not discussed here. Thus, like 
radiation, we have two types of pressureless matter (i) baryonic matter, 
which is formed due to nucleosynthesis and baryosynthesis of elemen- 
tary particles produced due to the decay of curvature inspired infiaton 
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and (ii) curvature-induced non-baryonic matter identified as dark mat- 
ter. After sufficient expansion of tlie universe, tlie matter dominates 
over radiation causing decelerated expansion as ~ t^^^ in tlie late uni- 
verse. 

It is interesting to note that, in the late universe, curvature-induced 
phantom terms in the Friedmann equation are Pph[l — Pph/2A] with pph 
being the phantom energy density. Moreover, curvature causes another 
constant A, which is analogous to negative brane-tension in Randall- 
Sundrum II theory of brane-gravity [15, 16, 17]. As a remark, it is nice 
to mention that these type of terms also appear in Preidmann equation 
based on the loop quantum gravity [18]. As branc-theory prescriptions 
are not used here, the curvature-induced A appearing in this model is 
identified as cosmic tension like the Refs. [13, 14]. 

Further, it is found that due to continuous increase in phantom 
energy density in the expanding universe, dominance of curvature- 
induced phantom begins at the red-shift z = 0.303 causing a cosmic 
jerk. As a consequence, a transition from deceleration to acceleration 
takes place in the very late universe. It is found that the phenomenon 
of cosmic acceleration will continue in future too. 

In the dynamics of the future universe, the cosmic tension plays a 
very important role as it opposes expansion of the universe. This effect 
is significant on sufficient increase in phantom energy density. Due to 
growth in phantom energy density in the expanding universe, it is found 
that matter will re-dominate cosmic dynamics, when phantom energy 
density will grow to twice of the brane-tension. As a consequence, 
universe will decelerate. 
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Increase in pph will still continue with the increasing scale factor 
a{t). So, phantom terms in the Friedmann equation will be negative 
when p > 2X. It is found that expansion will stop on the growth 
of the scale factor upto a maximum value = 1.88 x 10^° at time 
tm — 3.45 X lO^^^o (^0 being the present age of the universe) and the 
universe will bounce causing contraction in the universe. As a result, 
matter energy density will increase rapidly such that energy density 
and pressure density diverge as well as a = at icoi — 3.62 x lO^^io- It 
means that big-crunch singularity will be caused and the universe will 
collapse. This result is based on the classical mechanics. 

It is imperative to note that energy density , pressure and curvature 
will be extremely high near the time of the cosmic collapse tcoi- So, large 
structures of the universe will get smashed to elementary particles. 
This state is analogous to the early universe. So, quantum gravity 
effects can not be ignored near tco\ as quantum effects are important in 
the early universe. As a result, quantum particles will get produced due 
to rapid change in topology of the space-time. It is found that cosmic 
collapse can be avoided due to the back-reaction of these particles and 
the universe will expand exponentially after t > tco\- Thus, this model 
predicts revival of the state of early universe in future when t > tcoi- 

The paper is organized as follows. In section 2,/(i?)-gravitational 
equations are derived and trace of these equations are obtained. Curvature- 
inspired Friedmann equation is obtained in the early and late universe 
are obtained from the trace of gravitational equations. The early uni- 
verse is discussed in section 3. In this section, it is found that curvature- 
induced dark energy mimics quintessence and the early universe under- 
goes the power-law inflation. During curvature-induced quintessence 
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scalar behaves as an inflaton. Thus, the source of inflaton is known, 
which is an important consequence of this modeL Further, it is dis- 
cussed that the curvature-inspired inflaton decays during this period. 
This decay causes production of elementary particles and radiation, 
which thermalizes the early universe upto a very high temperature 
~ lO^^GeV. In section 4, it is discussed that the standard cosmology 
is recovered after exit from inflation and decelerates being driven by 
radiation and subsequently by matter. This section also contains tran- 
sition from deceleration to acceleration driven by curvature-induced 
phantom. Re-dominance of matter in future, cosmic collapse in future 
and its possible avoidance using the quantum gravity are addressed in 
section 5. In the last section, salient features of the present model and 
concluding remarks are given. 

Natural units(A;B — h — c — 1) (where ks, h, c have their usual 
meaning) are used here. GeV is used as a fundamental unit and we 



where G = Mp'^{Mp = lO^^GcV is the Planck mass), a is a dimension- 
less coupling constant, /3 is a constant having dimension (mass)^"^'') (as 
R has mass dimension 2) with r being a positive real number. 

Using the condition SS/Sg^" — 0, the action (2.1) yields gravitational 
field equations 



have IGeV"^ = 6.58 x IQ-^^sec and IGeV = 1.16 x 10^^ K. 



2. f{R)- gravity and Friedmann equations 



The action is taken as 




(2.1) 
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-/3(2 + r)i?(^+'-)i?^, = 0, (2.2) 

where stands for the covariant derivative. 

Taking trace of (2.2), it is obtained that 
R 



with 

□ = 



- &aUR - 3/3(2 + r) + /3ri?(2+^) = (2.3) 
1 d / , „„ d 



(2.4) 



In (2.3) 

□^(i+r) _ (i + r)[i?''ni? + ri?("-^)v^i?V^i?]. (2.5) 
Prom (2.3) and (2.5) 
— -^-[6a + 3/3(l + r)(2 + r)i?nni?-3/3r(l+r)(2+r)i?("-^)v'^i?V^i? 

+/3ri?(2+'") = (2.6) 

In (2.6), [6a + 3/3(1 + r) (2 + r)R^] emerges as a coefficient of UR due 
to presence of terms aR? and jSR^"^^^^ in the action (2.1). If a = 0, 
effect of R^ vanishes and effect of R^"^^^^ is switched off for (3 — So, 
an effective scalar curvature R is defined as [11] 

7^'- = [6a + 3/3(1 + r) (2 + r)i?n, (2.7) 

where 7 is a constant having dimension (mass)"^'' being used for di- 
mensional correction. 

Using (2.7) in (2.6), we have 

^J_yiA _ (7/r)^'-r(V'-2)[rnr + (i/r - i)v'^rv^y] 

-3(/3/r)(l + r)(2 + r)y('/'^-')v'^yv^y + /3ry('+'')/'' = 0, (2.8) 
where 

iRT - 6a 

Y^R'^ J J (2.9a) 
3/3(1 + r)(2 + r) ^ ^ 
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(2.8) is simplified as 

^ y- (7/r)^''[ny + (1/r- i)y-Vyv^y] 



-3(/?/r)(l + r)(2 + r)V^FV^F + /5rr(^/''+^^ = 0. 
Using (2.9a)in (2.9b), it is obtained that 



(2.96) 



6a 



+ [2,0'r{l + r){2 + r)h^]k 



[6q; - 7i?' 



(2.10) is re-written as 



-(1-r) 



Go; 



3/3(1 + r) (2 + r) 
+ Di? + (r - 1)R-^\I^'R\I^R 



= 0. (2.10) 



^ ^V^RV^^R + rR-^V^'RV^R 

Qa — ^R^ 



+r)(2 + r)/7']i?' 



,21 D2r-1 



7i?'' - 6a 



(l/r+2) 



= 0. 



(2.11) 



-3/3(1 + r)(2 + r). 
Experimental evidences [22] support spatially homogeneous and flat 

model of the universe 

dS'' = de - a\t)[dx'' + dy"" + dz""] (2.12) 

with a{t) being the scale factor. 

For a(i), being the power-law function of cosmic time, R ~ a~"^. For 
example, R ~ for matter-dominated model. So, there is no harm 
in taking 

^=4' (2-13) 
a" 

where n > is a real number and ^4 is a constant with mass dimension 
2. 

Connecting (2.11) and (2.13), it is obtained that 
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/5 



-1/3 



taking (— = — and ignoring complex roots as these roots 
lead to unphysical situations. Now, we have following two cases. 

Case 1 : The Ecirly Universe 

In this case, a{t) is very small, so (2.14) is approximated as 



a 

- + 
a 



2 — n — nr 



-1/3 



n(7A'-a-"^)2[3r(l + r)(2 + r)Y+yr 
x[6Q;-7^''a-"n'^'^' (2-15) 



as 



6a — ^A^a 
Integration of (2.16) leads to 

2 B 



-1. 



(2.16) 



2/3 



\aJ a(2+2A^) n(7A'-)2[3r(l + r)(2 + r)]i+V'-a(2+2A^) 

X y"a(^+'^+'"'-)[6a-7^'-a-n'+'/'' (2.17) 



with 



M ^ 2 - n - nr. 



(2.18) 



Case 2 : The Late Universe 

In this case, a{t) is large, so (2.14) is approximated as 

6aa^ 



+ 



2 — n — n{r — 1) — nr 



- 1 



n(7A^a-"0^[3r(l + r)(2 + r)]i+V' 



1677^7^'' L 

(6a)2+V^[a2"'- - (2 + l/r)7A"a"T 

(2.19a) 



as 







6a — ^A^a" 

for large scale factor a. So, (2.19a) is re- written as 

2 



+ 



(2.196) 
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where 



(2.20a) 



and 



/ 6q; \2 



1 



[3r(l + r)(2 + r)]-^-^/V 6 



n 



(2.19b) is integrated to 



(-) =- 

Va/ a 
with 



B 



2D 



(2+2JV) [2 + 2N + nr) 



^_ E{2 + 2N + nr) ^.^ 
~L)(2 + 27V + 2nr)"' 



(2.21) 



= 2 - 2nr. (2.22) 

Further, it is found that if M = 1, the first term on r.h.s. (right hand 
side) of (2.17)gives radiation. Moreover, if = 1/2 the first term of 
r.h.s. of (2.21) gives matter. So, setting M = 1 in (2.18) and N = 1/2 
in (2.22) to get aviable model, it is obtained that 

(2.23) 

1 



3 

nr = -, 



and 



n = 



r = 3. 



(2.24) 



(2.25) 



3. Power-law inflation, pcirticle creation and thermalization 

in the eEirly universe 

3 (a). Power- law inflation 

The approximated Priedmann equation (2.17), in the case of the 

early universe, looks like 
'd\2 B 8/3-V3 



/ay _ ±} 
\a) 



(7^3)2[l80]4/3 



7/l»a-Vi]'/3ia (3.1) 
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using definitions of M and as well as (2.24) and (2.25). In (3.1), 



/ 



- J. J. J 



X 



j a''/'{6a - -fA'a-'/^y/'da. (3.2a) 
It is noted that for 

/ \ 4/3 

a < I 'f Ay 6a] = (3.26) 

terms within bracket and the integral on the right hand side of (3.2a) 
are of the order of a^^/^. 
So, 



/ 



Thus, using (3.2a,b,c),(3.1) is approximated as 

2 B 16f3-^/^ 



0' 



ll(7A'^)-V3[i80]4/3 



3/2 



a-3/4 _ ^-3/4 



7/3 



(3.2c) 



(3.3) 



It is interesting to see that a radiation density term B/a* emerges 
spontaneously. This type of term emerged first in brane-gravity in- 
spired Priedmann equation. So, analogous to brane-gravity, here also 
B/a^ is identified as dark radiation. Other terms on r.h.s. of (3.3) are 
caused by linear as well as non-linear terms of curvature in the action 
(2.1). These terms also constitute energy density term 



in 



G 



16/3- 



-1/3 



,3/2 



-3/4 



6a 1 7/3 



(3.4) 



.ll(7A3)-V3[180]4/3 

(taking real root of (— 1)"^/^ as above) satisfying the conservation equa- 
tion 



Pde + 3-(pde+Pde) = 0. 



(3.5) 



Connecting (3.4) and (3.5), equation of state (EOS) is obtained as 



(3.6a) 



where 
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/ 



(3.66) 



(3.6a) is the scale factor-dependent equation of state parameter ,vahd 
for ap < a{t) < Uc- Such an equation of state parameter is obtained in 
[10] also. It yields 

pZ+pZ>o and pI: + SpZ<o, 

for ap < a{t) < ac- It shows that DE, having energy density (3.4) 
mimics quintessence dark energy [9, 10, 11, 12, 13]. 

Here investigations start at the Planck scale, where DE density is 
obtained around lO^^GeV^. So, (3.4) is re-written as 

pll = Fa3/2[a-3/4 - af/^y/^ (3.7a) 

with 



75^-3/2 



lO'^a 



Thus , (3.6b) and (3.7b) imply 



-3/4 _ 3/4 



-7/3 



f = F. 



(3.3) and (3.7a) imply 



d\2 B 



i4 ■ 3 
using G = Mp^ = 10-38GcV"l 



StT X 10^7 / a \ 3/2 pQ-3/4 _ ^-3/4^ 7/3 

ap/ 



-3/4 



-3/4 



(3.76) 



(3.7c) 



(3.8) 



As Op is expected to be extremely small, so dominates over the 
radiation term in (3.8). Moreover, (3.7a) shows that vanishes at 
a — ac- So, for Op < a{t) < ac, cosmic dynamics is given by 



d\2 _ 87r X 10^^/ a \3/2 
aJ 3 Vop/ 

_ 87r X 10=^^ / a \ -V4 
3 Vop/ 



,-3/4,^-3/4 7/3 
-3/4 -3/4 
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(3.9) 



(3.9) integrates to 



a{t) =ap 1 + -^^{t - tp) (3.10) 
8v37r 



showing acceleration as a > 0. 

If expansion (3.10) yields sufficient inflation in the early universe, 

^ = 10^^ (3.11) 

ap 

The universe comes out of the inflationary phase at t — tc when a{t) 
acquires the value Oc. So, from (3.10) and (3.11), it is obtained that 

~ 7.77 X lO^tp (3.12) 

using (3.11). 

3(b). Realization of curvature induced quintessence dar\s. 

energy through scalar field 



In 3(a), it is found that curvature induced quintessence dark energy 
causes power-law inflation. Though, we have a gravitational origin of 
quintessence DE here, it is natural to probe a scalar 0(i, x) giving pde 
obtained above using the scheme adapted in [12, 19]. With V{(j)) as po- 
tential and having minimal coupling with gravity, in the homogeneous 
and flat model of the universe (2.12), 0(t,x) obeys the equation 

□0(t,x) + y'(0) = (3.13a) 

and has mass dimension equal to 1. Here, T^'(0) = dV/dcf). 

In the homogeneous and flat model of the universe (2.12), (3.13a) is 
obtained as 

0o + 3-0o + l^'(0)^=0o = O (3.136) 
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as x) = (po(t) due to homogeneity. DE density and pressure for 
4>{t) are given as 

Pde = + Vi<t>o) (3.14a) 

and 

Pde = - Vi^o) (3.146) 

It is interesting to see that conservation equation (3.5) yields (3.13b) 
for DE density and pressure of the (po{t)- fluid given by (3.14a) and 
(3.14b), 

It is demonstrated above that quintessence DE, given by (3.7a) and 
(3.7b), derives the power-law inflation with a{t) given by (3.10). The 
Raychoudhuri equation (which is obtained connecting Priedmann equa- 
tion and conservation equation) yields 



H = -47rG'(pde +Pde) = -^nG<j)l 



(3.15) 



using (3.14a) and (3.14b). H, in (3.15), is obtained from (3.10) as 



„ a Mp 
H — - — — 
a V37r 



1 + -^(t-tp) 



8V37r 



(3.15) and (3.16) yield 

00 



Ml 



47rv^ 



i + -^(t-tp) 



8V37r 



(3.17) integrates to 



00 



Mp\\-ln 

TT 



8V37r 



V 327r \ap) 



(3.16) 



(3.17) 



(3.18) 



using (3.10) and 0^ = at Planck scale. Further, connecting (3.11) 
and (3.18), it is evaluated that 



0^ = 6.43Mp 



(3.19) 
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being 0o at the end of inflation. (3.18) yields the relationship 

a{t) = apel"^"^^'^]. (3.20) 

As 00 is the quintessence scalar, giving curvature induced quintes- 
sence DE. So, using (3.6a,b) and (3.7a,b,c), it is obtained that 



^(0o) = ^(pde-Pde) 



5 



4 



3/2^.-3/4 -3/4 X 



7 



30 



(3.21) 



with a{t) given by (3.20). As V^(0o) = at 0o = (t>o- it shows that quin- 
tessence scalar 4>q falls from the high hill of the potential to the ground 
state at the end of inflation. For 0o < 0O) (3-21) is approximated as 

V{(l>o) ~ ^Fa-'/'e-[<^°^^'^] (3.22) 

with F given by (3.7b). Thus, it is obtained that curvature induced 
quintessence DE driving power-law inflation can be realized through 
quintessence scalar (f)o{t) given by (3.20). 

In non-gravitational models of DE, origin of quintessence is not 
known. This model has an advantage over non-gravitational DE models 
due to gravitational origin of quintessence scalar. 

3(c). Pcirticle creation during inflation 

Creation of spinless bosons 

4>o{t) is the background field deriving inflation and 4>{t,x) can be 
realized as 0(i,x) = 4>o{t) + 50(i,x) with 50(t, x) being the quantum 
fluctuation. Here, perturbations in the metric components are ignored 
for simplicity. So, from (3.13a) and (3.13b), we obtain 

nS(l){t, x) + +V"{<f))^=^,S(l){t, x) = (3.23a) 
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In the space-time (2.12), (3.23a) looks like 

Q2 Q2 Q2 



a i 



+ 



dx^ dy^ dz^- 



5(f)(t,x) + 3-5(j)(t,x) - a-^ 
a 

Using the decomposition 



+ 



(3.236) 



50(t,x) = Yl[Mt)ake''-' + <l>l{t)aie-''-% (3.24) 

— oo 

a{t) given by (3.10) and V{(f)) given by (3.22), (3.23b) is obtained as 

72a. OA ^A. ri92vrA;2 60i 

(3.25a) 



dVfe ^ 24 mpk 



drf 



for a < Op. Here 



50(?7, x) and 11 



r] dr] 
1 + 



Mp 



+ — 







{t - tp) 



(3.256) 



8V37r 

,x^)dri in (3.24) satisfy quantum conditions 



[5(/)(?7,x),50(?7,xO] = O 

[n(r/,x),n(r/,x')] =0 

[50(r?,x),n(r;,x')] = z(-5)-V^5^(x-x') 
with 5^(x — x') being the Dirac delta function. 

(3.25a) integrates to 

(l>k{t) = (pkiv) 



(3.25c,d,e) 



-23/2 



(3.26a) 



CiJir{-kVl927rr]-y7apMp) 
+C2Yir{-kVl92nr]-y7apMp) 
where Jn{x) and Yn{x) — (— l)"J„(x) are Bessel's functions with 

(3.266) 



n — ±ir — ±-. 



60 



11.52 



with i = ^/—l. Ci and C2 are integration constants. Moreover, from 
(3.12) and (3.25b), it is obtained that 1 < rj < 7.77 x lOHp. So, for 
ap 10~^^, which is possible for ap being the scale factor at Planck 



scale, \k\/1927rr] '^/7apMp\ is large. 
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Using identities 



dJ-n{x) dJn{x) 2sin{n7r) /o o« ^ 

JnW — 1 1 J-n{x) = , (3.26c) 

ax ax TTX 



Y^{x)^-^^ - ^Y_^{x) = _^£M!^, (3.26d) 
and (A6), 

= CiJir(-kVl92^r]-y7apMp) 

and 



are normalized. As a result, Ci and C2 are evaluated as 
For large x, 



(3.26e,/) 



and 



7apMpr]lsin{r7r) . 
JJx) ~ 4/ — cos[x - 7r/4 - n7r/2] (3.27a) 

V TTX 



/ 2 

y;,(,T) ~ W — sin[x - T^/A - 7171/2] (3.276) 



Using approximations (3.27a) and (3.27b) in (3.26a), we obtain 



lapMp _g 



0L.(?7) ~ 4 / i^^Tj — iCiCosix — Ti / A — im 12) 



+iC2sin{{x -7r/4-ir7r/2) , (3.28a) 

where 

^ ZM^. (3.286) 

cosx and smx have similar characteristics as both are periodic func- 
tions. So, there is no harm in taking arbitrary constants Ci and C2 as 
Ci = iC2- With this identification, (3.28a) is obtained as 



47rv37r/c 
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with 

ra.M.,?nV^_ (3.296) 

It is obtained using the normahzation condition (0^, 0fe') = ^kk' with 
scalar product (A6) at the hyper surface r] — etai. 

Here, in— and out— states are obtained when fj tends to — oo and 
oo respectively. So, 

oo 



— oo 
oo 



= 5^[0r(^)«re^'-^>0r(^)«fe-^^-1 (3.30) 

—00 

0fe"*(^) cind (p'^'iv) s-rs related through through Boglubov transfor- 
mations (A7), where Boglubov coefficients ak and (5k are defined by 
(A8) and (A9) respectively and satisfy the condition (AlO). 

Using the above definition of in— and out— states, we have 



47r\/37rfc 



and 



</,-(r^) = C2W '^'''^j-f' e-'-"/^r^-^e'("^^'^"'/^°^^^+"/^) (3.316) 
47rv37rfc 



Using (3.31a) and (3.31b) in (A8) and (A9), Boglubov coefficients 
are obtained as 

ak = -4iC|77j-^^[8r7r^ - TiQle-^^^^^i"" (3.32a) 



and 



where 



and 



(5k = 32iC|e^'^/277ri^ (3.326) 
C3 = C2e-''"/V7apMp/47rA;7r (3.32c) 
C4 = 8kVSn/7apMp. (3.32d) 
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Introducing (3.32a) and (3.32b) in conditions (AlO), it is obtained 



that 



yielding 



28C|C4^r^^ = 1 (3.33a) 



C|e-''^ 11/16 

TT 



(3.336) 



So, 



with 



\ak\^ = l + IO24C3V'' = 1 + ^ (3-34a) 



A = 1024[C2e-'''^/V7apMp/47r7r]^r7f^^ (3.346) 

Further non- vanishing obtained from (3.32b), shows creation 

of spinless quantum particles during inflation. 

Moreover, the rate of creation of scalar particle- antiparticle pairs per 

unit time per unit volume [24] , due to decay of 0, is obtained as 

00 00 
r^_.#$ = -ln{ H \ak\-^)/V^ = Yl ^^l«felV^4 

K=—oo K=—oo 

^ E = ic(2)M 

K=—oo 

^ 25MpiC(2) 
where Riemann zeta function ^(2) — and 



(3.35a) shows that creation rate of these particles decreases as time 
increases. 

Creation of spin-1/2 Dirac fermions 



If il>{t,y:) are spin-1/2 Dirac spinors (being mathematical represen- 
tation of spin-1/2 elementary particles) during this phase , it is natural 
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for these quantum fields to interact with curvature induced quintes- 
sence scalar with interaction term —h'ijj(j)oip{ijj = where i/j"^ is 
the adjoint of ip and 7° is the time-component of the Dirac matrix 
in flat space-time. Here h is the dimensionless coupling constant. In 
the present model,it is assumed that few spin-1/2 elementary particles 
exist during this period, but these are not sufficient to effect cosmic 
dynamics. Rather, these particles act as seeds for creation of more 
particles during the inflationary period. 
Prom the action 

= J d'^x^/-9i'[iYD^ - h(f)o]i', (3.36a) 

the Dirac equation is obtained as 

i-f'^Di.i; - hM = (3.366) 

where is defined in (All). 

In the space-time (2.12), (3.36b) is obtained as 

[a{t)fdo + Tda + ih(l)oa{t)]iP = (3.36c) 

with a = 1, 2, 3. Operating (3.36c) by [a{t)^^do + ^"'da — ih(l)oa{t)] from 
left and using the decompositions of tp given in Appendix A, we get 
the differential equation for fk,s{gk,s) as 



• . 72 /, 

/+ -/+ [- + ^e-(0oa + a0o) + /iVo]/ = 0, (3.36d) 

a a^ a 



where / = fk,s{9k,s)- 

(3.10) shows that for t < 25.56tp, 



T ^ ^^{t - tp) < 1. (3.37a) 

v37r 



So, for t < 25.56tp, a{t) (given by (3.10)) is approximated to 

a{t) = a(T) ~ ape''. {3.37b) 
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(3.20) and (3.37b) imply 

for t < 25.56tp. 

Using (3.37a), (3.37b) and (3.37c), (3.36d) is obtained as 



(3.37c) 



d'f , df 



+ ^ + 



4M| 




;i - 2r + 2t^) + ieh\l^{l + r) + ^h\^ 



32 



32 



(3.38) is obtained in the simplified form as 



d'^f , df 



+ -1^ + 



df^ df 



M2 + 



2f 



/ = 0, 



where 



alMl 



247rA;2a|,M|, + 3h'^alMj,/8' 



(-67rfc2 + ie/ia|M2 ^3/32) 
^Q'Kk'^alMl + 3/i2a|>iV/|,/32 



and 



with 



f = Mr 



M = L-2[-67rA;2 + ie/ia|M|, ^3/32]. 
(3.39a) yields the solution 

Using (3.40) in (A14) and (A15), it is obtained that 
with Us and {is given in Appendix A. 



/ = 0, 
(3.38) 

(3.39a) 

(3.396) 
(3.39c) 
(3.39ci) 
(3.39e) 

(3.40) 

(3.41a) 
(3.416) 
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Using in- and out-statcs defined as above, we obtain ipfl ^ and V'/fc's 
for f < and f > respectively. Similarly, _g-^ and ipjY{-k -s) 

are obtained. 

In this case, the Bogolubov coefficients ak,s and Pk,s, defined by (A17) 
and (A18), are obtained as 



(3.42a, b) 



as these satisfy the condition (A16). 

When t > 25.56tp, a{t) (given by (3.10)), ^{t - tp) > 1 and a{t) 
is approximated to 

Mp 



a{t) ~ ap 



as well as (3.20) and (3.43a) imply 

Mf 



(t - tpf. 



(3.43a) 



00 ~ 8^^[ln(t -tp) + ln{Mp/VS7r)] 
V 327r 



Mp 

'V32^ 
So, (3.36d) looks hke 



8-^[-{t -tp) + l + ln{Mp/V3Ti)] (3.43b) 



/ + 



8 



+8ieh 



(t - tp) 
Mp 





r Vs^ 1 


-ap 


[Mp(t-tp)\ 



VS27r{t - tp) 



f + 



(9 + 8/n(Mp/V37r) - 8{t - tp)) 



+2h'—^{[{t -tp)-l- ln{Mp/V37i)y 



/ = o. 



(3.44a) 



Now, for 25.56tp < t < 7.77 x lOHp, (3.44) is approximated to 



/ + 



{t - tp) 



f + 





r \/3^ 1 


-Op 


lMp{t-tp)\ 



16 Ml 1 ~ 

+ 2h^—^{t-tpf f = 0, (3.446) 

TT J 



For k small, (3.44b) is approximated as 
/ + 

which integrates to 



f^C{t- tp)-^/^J^r/4{±hMp{t - tpf/\Ph:). 



(3.44c) 



(3.44d) 
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For k large, (3.44b) is approximated as 



k 



^371 



16 , 



which integrates to 

/ = dt - tp)-^/2 J_i/2(-76(t - tp)-'). (3.44/) 

Solutions (3.44d) and (3.44f ) also yield Bogolubov coefficients ak,s and 
f3k,s given by (3.42a,b). 

Results (3.42a,b) show that elementary spin-1/2 particle-antiparticle 
pairs are created during the inflationary prod. The rate of the creation 
is obtained as 

oo oo 
s=±l fc=— oo s=±lfc=— oo 

n 2 



= 4Zn2 /V4 = 100Mpln2 /a^VSn l + Mp{t-tp) /sqrtSn 5 (3.45) 

using C(0) = —1/2 obtained through analytic continuation. 

Fluctuation of <j)o{t) during inflation and thermalization 

(3.13b) is an equation for damped harmonic oscillator (f)o{t), which is 
the curvature-induced quintessence scalar given by (3.20). It indicates 
small oscillation around (poit) at every instant of time during its fall 
from high hill of the potential to the ground state given by ^(^c) — 
0. For example, at t — t, the background scalar is 0o(^)- So, small 
homogeneous fluctuation S(f)o{t) (it is taken homogeneous being non- 
quantum and due to homogeneous space-time ) around (f)o{t) — (f)o{i) 
will satisfy the equation 



SMt) + 3-50o(t) + V^"(0)^=^„(,-)50o(t) = (3.46) 
which is obtained introducing 0o(^) — 0o(^) + ^<Po{'t) in (3.13b). 
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Further, incorporating (3.10) and (3.22), (3.46) is obtained as 

<i^^24<H^^ (3.47a) 

df)^ T) df] 

with 

= 2407r X io37e^o(t)M-i^_ ^3 47^^, 

(3.47a) integrates to 



SMt) = C^"'^/'^-23/2(&ry) ~ y/2/^7]-^^cos{bri - tt/2 + 237r/8). 

(3.48) 

(3.48) shows that fluctuation S(f)o{t) oscillates around 0o — 4>o(t) 
with amplitude decreasing with time. So, energy will be released as 
radiation due to these fluctuation (as (amplitude)^ yields energy for 
the classical held). 

As this phenomena will occur at every instant of time during the 
inflationary period tp < t < tc, the continuous release of energy during 
fall of 00 will thermalize the created elementary particles. The energy 
released, during the inflationary period, has density 

V{0) - T/(0e) ~ lO^^GeV^ (3.49) 

It shows that the released energy will thermalize created elementary 
particles upto sufficiently high temperature Tc such that 

15 " 

It gives 

Tc = 4.8 X lO^^GeV. (3.50) 

Like [12], here also, this radiation is recognized as cosmic microwave 
background radiation (CMB). 

Thus, from the above analysis, it is obtained that, during the infla- 
tionary period, curvature-induced quintessence scalar causes creation 
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of elementary particles and energy released in the form of radiation 
heats up these created particles upto very high temperature. So, these 
particles are highly relativistic having thermal equilibrium with the 
emitted radiation. 

4. Deceleration driven by radiation and matter as well as 



Deceleration driven by radiation 

(3.4) shows that curvature-induced quintessence energy density van- 
ishes at a = ttc. As a consequence, universe exits from the inflationary 
phase. During this period, dark energy causes elementary particle - 
antiparticles as well as CMB thermalizing the created particles upto 
temperature = 4.8 x lO^^GeV given by (3.50). 

Moreover, in (3.3), we get curvature-induced dark radiation term. 
Thus, we have two sources of radiation. As temperature of CMB, ob- 
tained here, is very high, dark radiation too will have thermal equilib- 
rium with CMB. So, energy density of created particles, dark radiation 
and CMB together will have energy density 



Thus, at the end of inflation (a = Cc) , we recover standard model of 
cosmology and (3.3) reduces to 



Acceleration in late and future universe 




(4.1) 




(4.2) 



(4.2) integrates to 



a 



{t) = ac[l + 4Mp^yn/15at{t - Q]'^^. 



(4.3) 



Matter-dominance and deceleration 
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In the late universe, the effective Friedmann equation is given by 
(2.21). Using (2.23)-(2.25) in (2.21), we obtain 



ay 
a) 



C 8D 



+ 



,3/4 



15 



5^ 
6d' 



,3/4 



where 



and 



-3/4 



D = 



A-kG 

-3/2 

ac 



287rG / 1 
135 V180/? 



1/3 



-1/4 



135 V 180/3/ " . 



(4.4) 



(4.5a) 



(4.56) 



being obtained from (2.20a) and (2.20b) using (2.23)-(2.25) and (3.2b). 

The first term, on r.h.s. of (4.4), emerges spontaneously and has the 
form of matter density, so it is recognized as dark matter density like 
dark radiation. 

Before proceeding further, it is useful to remark that, like radia- 
tion, we have two types of matter too (i) dark matter given as ^^Q^i 
in (4.4a), which is non-baryonic and (ii) baryonic matter formed by 
elementary particles ( produced during inflation) through various pro- 
cesses of standard cosmology such as nucleosynthesis, baryosynthesis 
and recombination of hydrogen (not being discussed here) . 

According to WMAP results [23], present density of pressureless mat- 
ter (baryonic and non-baryonic) is obtained to be p"^^ — 0.27p^ and 
present dark energy density p^^Q = 0.73pq^ with 

^0^ 



SttG" 



where current Bubble's rate of expansion Hq — lOOhkm/ Mpcsecond — 
2.32 X 10-^2/iGeV and h = 0.68. Thus, 

2.9 X 10-^^GeV^ 



Po 



Ho = 1.58 X 10-^^/iGeV = [0.96io]~^ 



(4.10a) 
(4.106) 
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and 

to = 13.7Gyr = 6.6 x lO^^GeV-^ (4.10c) 
Using these values, it is obtained that 

(mat) ^ 3C 0.27H^ 

with the present scale factor normalizing as 

ao = 1 (4.116) 

and 

p£ = ^a^'/^ = 0.73i/oV/^ (4.12) 

from (4.6). 

Connecting (4.4), (4.7) and (4.8), it is obtained that 



^ + 0.73i^„V/^{l-°-y} (4.13) 
(4.13) shows that 



d\2 _ 0.27//n 
aJ 



for a < 0.767 and 
for a > 0.767. 

It means that a transition from matter- dominance to DE-dominance 
takes place at 

a* = 0.767 (4.14) 

giving red-shift 

= — - 1 = 0.303 (4.15) 
a* 

which is very closed to lower limit of given by 16 Type supernova 
observations [2]. Thus, for a < 0.767, (4.13) is approximated as 

(ty^ojif (4.16) 
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using matter density given by (4.11a) and 

[0.73p^'-a3/4]2 « 0.73p^'a^/^ 

for a being less than 1. 
(4.16) integrates to 

a{t) = ad[l + ^Vo:27Hoa-'^\t - ta)f^ (4.17) 

(4.13) is approximated as (4.16) when 

» 0.73Hy/\ 

but as a{t) approaches very close to a = a*, 

^-^^ ^ 0.73Hy/' 
a"* 

So, (4.16) needs to be modified as 
which integrates to 

3 
2 



(4.18a) 



a{t) = ad[l + ^y/OMHoa/^\t - td)^/^ (4.186) 



It shows decelerated expansion as a < 0. In (4.17), td —— 386kyr = 
2.8 X 10~Ho (WMAP result) is the decouphng time and the scale factor 

at t = td is given by l/ad = 1 + Zd^ 1090 with Zd = 1089 (WMAP 
result). Using these values and a* from (4.14), (4.10b) and (4.17) yields 

U = 0.63^0- (4.18c) 



Phantom dominance and late acceleration 



When a > 0.735, (4.13) is approximated as 



-3/4 



2A 



(4.19) 



with Hq given by (4.10b). 
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(4.19) integrates to 
a(t) = 



2A 



2A 

2-1 -4/3 



--«„Va73((-i,)}] (^^g^j 

as a*^^^ = 1.22. (4.20a) shows and acceleration and it is singularity- 



free. 

Prom (4.20a), it is obtained that 



1.7p% 11 



-3/4 



(4.206) 



This shows a > , when 



. 11-3 



A ' 3 

Thus, solution (4.20a) yields acceleration in the late universe. 
(4.19) shows that accelerated expansion (4.20a) stops at a = Og sat- 
isfying the condition 

0.73p°,a^/^ = 2 A. 



(4.21) 



Further, (4.20a) yields 
1 = Oo = 



2A 



+ 



0.73p; 



2A 

2-1 -4/3 



Using (4.20a) and (4.10b) in (4.22), A is evaluated as 

A = 1.17p,V 

a{t), given by (4.20a), acquires the value by the time 



(4.22) 



(4.23) 



ie = t, + -[i7oVa73]-Yl.22 
= 0.63io + 2.49io = 3.12^0 



0.73p| 



2A 



(4.24) 
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5. Re-dominance of matter, collapse in the future 
universe and its avoidance 

(a) Re-dominance of matter and deceleration in future universe 

It is shown above that phantom-dominance ends at t — te- As a 
consequence, matter re-dominates and (4.13) reduces to 



a3 



(5.1) 



(5.1) integrates to 



a{t) = ae[l + ^]Vo:27Hoaf/^{t - te)f^ (5.2) 

showing decelerated expansion as a < 0. 
Moreover, though at a = Og 

vanishes, it will be negative for a > a^- As a result, at a certain value 
of a{t) being > a^, a — 0. Thus am is the maximum of a{t). It is 
because, for a > a^, [a/a]^ < leading to imaginary expansion rate 
a/ a. So, ttm satisfies the condition 

„.,.|0^_^|^27 (5.3„) 

(5.3a) yields 

am = 0.86/[p^'"]2/9 = 1.88 x 10^° (5.36) 
with Pq given by (4.10a) 

Connecting (4.21), (4.23), (4.24), (5.3a) and (5.3b), it is obtained 
that 

t„ -3,45x10%. (5.3c) 
(b) Contraction in the future universe and its collapse 
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Further, it is interesting to note that the curve a = a{t) will be 
continuous att = tm, but the direction of tangent to this curve (pointed 
ai a — Qm) will change yielding d < for i > i^. It means that 
universe will retrace back at t — tm and will begin to contract for 
t > tm- During the contraction phase, a{t) will decrease with time. So, 
term proportional to a^^ will dominate over terms proportional to a^^^ 
in (4.13). As a result, we have 

H^^ ^0.27H^a-\ (5.4a) 

yielding 

= (^) ~ -VMlHoa-'/^ (5.46) 

(5.4b) integrates to 

ait) = am[l - ^VMlHoa-^/'it - t^)]'/' (5.5) 

showing decelerated contraction as d < 0. 
(5.5) yields a{t) = at 

i — i — + _i_ ^ ff-in^/s 

t — tcol — "T 27 

= 3.62x10^%. ^gg^ 

So, at t — tcoi, dominating energy density term 

p{mat) ^ o.27p^7a3, (5.7) 

in (5.1), will be infinite. These results show cosmic collapse at this 
particular epoch. 

In what follows, a possibility to avoid the cosmic collapse, obtained 
by the classical mechanics, is demonstrated. 

(c) Scaleir field for matter and peirticle creation in the future 
universe 
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In section 3(b), curvature-induced quintessence dark energy is real- 
ized through the scalar 0. Using the same approach, matter term (5.7), 
in Friedmann equation (5.1), can be realized through a different scalar 
$(x, t) obeying the equation 

□ $(t,x) + 1/'($) =0. (5.8a) 

Here, V'{^) — dV/d^ and ^{x,t) has mass dimension equal to 1. 

In the homogeneous and flat model of the universe (2.12), (5.8a) is 
obtained as 

'^{t) + 3-6(t) + l^'(*)$=$(t) = (5.86) 
a 

as $(i,x) = $(t) due to homogeneity. 

Density and pressure for $(i) are given as 

and 

= = _ y^^^) ^5 9^^) 

Connecting (5.4), (5.7), (5.9a) and (5.9b). it is obtained that 

- ~ -V87r/3Mpi$. (5.10a) 

(5.10a) integrates to 



a{t) = a^e-VS^^^'^*-*-). (5.106) 
(5.7), (5.9a) and (5.9b) yield 



cr 

Vi^) = 0.135^ 



cr 

3^35^0 ^V2^M-^{^-^^) 



(5.10c) 

It is discussed above that near t — tcoi, energy density will be ex- 
tremely high. This situation is analogous to the state of early universe 
with high energy density and large curvature. So, near collapse time. 
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classical mechanics is not the appropriate theoretical machinery and 
investigations should be done using quantum field theory. These ar- 
guments prompt us to resort to creation of quantum particles due to 
gravitational changes and probe its back-reaction on the future uni- 
verse. 

If x) is the quantum fiuctuation in homogeneous background 
field we have the equation of motion for 5^{t, x) as 

05^t, x) + |$=$(t)5$(t, x) = (5.11) 

obtained from (5.8a). 

In the space-time (2. 12), (5. 11) is re- written as 

(I 0'^ 0*^ 

x) + 3-J^t, x) -a-'\— + — + —1 5$(t, x) 



a 



dx'^ dy"^ dz^ 



+y"($)$=$(t)5$(i,x) = 0. (5.12) 
Using the decomposition 

oo 

5$(t, x) = Y,[^kit)ak^^-^ + $:(t)4e-'^^-^1 (5.13) 

— oo 

as well as results (5.10c) and a{t) given by (5.5), (5.12) is obtained as 

m) + ^{{-^Lr^'i + -^}*.«) = (5.14a) 

with 



fl=[l- -V^Hoa-JI\t - U]- (5.146) 

Here prime (') signifies derivative with respect to fj. 
(5.14a) integrates to 

Mv) = r^'^[CiJin{hfi'/^) + Yin{hfi'l% (5.15a) 

where 

h = JJ^^ (5.156) 
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and 



n 



3V7/2. 



(5.15c) 



Using identities (3.26c), (3.26d) and (A6), Ci and C2 in (5.15a) are 
evaluated as 



Ci — C2 



-3/2 



1/2 



(5.15d,e) 



-3V0^Hosinh{mr)Vfif'^- 
where V is the 3-volume. 

Using approximations (3.27a) and (3.27b) for Bessel functions, (5.15a) 
is obtained as 

^kiv) ^ Cr^-'/^e-'^/^e*''*'''" (5.16a) 

with 



C 



iTidii, ■ C""'- 



1/2 



(5.166) 



- Hosinh{mr) Vfjl'^ - 

Like section 3(b), here also in— and out— states are defined such 

that 



(5.17a) 



and 



(5.176) 

Using (5.17a) and (5.17b) in (A9), Boglubov coefficient is ob- 
tained as 

(5k = 47ra^/2e=''^^/^e-^'^/277r73sm/i(37rV7/2). (5.18) 
Further, (5.18) yields 

|/5fe|2 = 167r2a^e3"^/2^rV9sm/i2(37rv^/2), (5.19) 



which is non- vanishing. It shows creation of spinless particle- antiparticle 
pairs. 
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Energy density for the quantum fluctuation 6^{t, x), satisfying (5.12) 

is obtained as 

1 . * 1 21bH^ 

^ a 



1 

4 



/ dS^t, x) \ / dS^* {t, x) \ / dS^t, x) \ / dS^* {t, x) \ 
2a(t) l\ / V ) ^ \ ) \ / 

^ / 9M>(f.x) \ / g^$*(t,x) \i (5.20) 
V dz J\ dz J .' 

Incorporating (5.13), (5.20) yields 

^ rl . . /I 91 ^ff^ 1.2 1.2 

— oo ' 

(5.21) 

Now, renormalized energy density of created spinless bosons is ob- 
tained as 

* I IT) I * * I Ctuf^ I * 

Pcreated = < in\p \in > - < m\p |^7^ > 

oo . 

= -J20,27CWH',[-ala-' + 9ala-'] 



-oo 

oo 



— oo 

-2C(0)0, 27CWH'o[Ula-' + 9«^« 



9 ■ m J ^5 22) 

-2C(-2)|/?,n^ + 0.73aia-«] 

= ^^/0:277r3aL/^(f ) (e-^.m/.3(3.V7/2)^r^'/^«-^ 

using (5.17a), (5.17b) and (5.19), C(0) = -1/2 (obtained through ana- 
lytical continuation) and C(~2) = 0. 

(d)Back - reaction of created particles and escape 
from the cosmic collapse 

It is natural to think that created spinless particles (obtained above) 
will effect cosmic dynamics. As a consequence, FE (5.4a) is modified 
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as 



Stt 



(5.23a) 

when energy density of created particles will be significant. Here, 



X 



The solution of (5.23a) can be taken as 

a = acoiexp[\D{tcoi - t)}\ + 7|£'(tcoi - t)\'^], 



(5.236) 



(5.24) 



where Ocoi is the scale- factor att — tcoi, D is a constant of mass dimen- 
sion and 7 is a dimensionlcss constant. 
(5.24) yields 

H^-D[l + 2-f\D{t,oi-t)\]. (5.25) 

(5.24), being solution of (5.23), will satisfy (5.23). So, 

0.81 



—M'pD'[l + 4^\D{t,oi-t)\+4Y\D{t,oi-t)\\ - ^ 

OTT OTT 

X [1 - 3|D(teol -t)\- 37|D(teol -t)f] 

+Xa-J[l - 5|L>(icoi -t)\- 57|L'(icoi - t)\% 
(5.26) yields 



(5.26) 



Stt 



0.81 

"StT 



Am|.D27 = -^^-^MlHh:^.-bXa~' 

ZTT OTT 

^M],D^^^ = [-37+- 

2 



27r 



'P-'-'0"col 

0.81 



col 

2 u2„-3 



87r 



25 

57 



(5.27a, 6, c) 

comparing constants terms on both sides as well as coefficients of 
\D{tcoi-t)\ and \D{t,oi-t)\\ 

From (5.27b,c), it is obtained that 

0.81 



25Xa-J 



-9- 



87r 



■MlHla-l 



(5.28) 
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This equation shows that as universe contracts matter density will de- 
crease and the effect of back-reaction of created particles will increase 
after a certain epoch t — to- 
Using (5.28) in (5.27a,b), it is obtained that 

(5.29a) and (5.29b) yield 

Also (5.29a) shows that acoi can not vanish, because D will be diver- 
gent. Thus, it is found that back-reaction of created particles will make 
energy density finite at t = tcoi- 

Planck scale is the fundamental scale. It suggests the largest energy 
mass scale as Planck mass Mp. At this scale, energy density is obtained 
M^/Stt^. So, we have 

connecting (5.27a) and (5.28). 
(5.31) yields 

acoi = 2.25 X 10"^^ (5.32) 

and 

D^^. (5.33) 



6. Salient features and concluding remarks 

Here, a cosmological picture is obtained from the gravitational action 
containing the linear Einstein term as well as non-linear terms and 
R^. This is a gravitational action with no other field except the scalar 
curvature R. The approach of this paper is different from the approach 
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in Refs.[5, 7, 8]. This approacli iias an advantage to have power to 
explain 

(i) power-law inflation in the early universe and graceful exit from this 
phase, 

(ii) creation of SM particles, 

(iii) recovery of the standard cosmology with the cosmic background 
radiation with extremely high initial temperature ~ lO^^GeV, 

(iv) deceleration of the universe driven by radiation emitted during the 
inflationary phase and particles in thermal equilibrium with radiation, 

(v) deceleration driven by curvature induced dark matter and baryonic 
matter caused by various processes like nucleosynthesis, baryosynthe- 
sis, hydrogen re- combination of elementary particles created during 
inflation, 

(vi) dominance of curvature-induced phantom at red-shift z — 0.303 
(which is consistent with observational results), 

(vii) transient acceleration driven by phantom in the very late universe, 

(viii) re-dominance of matter, contraction of the universe, collapse of 
the universe at time icoi — 3.62 x lO^^^o , 

(ix) avoidance of collapse due to creation of particles near the time tcoi 
as well as its back-reaction and 

(x) rebirth of the universe after tcoi- 

Thus model gives the revival of the state of the early universe at time 
icoi — 3.62 X lO^^io . Interestingly, above investigations show that the 
Ricci scalar R has dual role as a geometrical field as well as a physical 
as physical and geometrical terms in the above theory are obtained 
from the scalar curvature R noted in earlier works [25]. 

Appendix A 
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Boson case 

In the case of minimal coupling of $ to gravity and having interaction 
with (f) as — (l/2)g'0^<l>*$, we have the equation 

x) + [g(P\t) + m|]$(t, x) = (Al) 

with m$ being mass of $(i,x). 
$ is decomposed as 

oo 

*(i,x) = J2i'^k{t)a,e^'-' + mWke-'H (^2) 

— oo 

In (A2), akiajj are creation (annihilation) operators satisfying quan- 
tum conditions 

K, dk'] = = [4, 4,] 

[ak,al]^{27r)%k' {A3, A, 5) 

The scalar product [24] is defined as 

($1, $2) = -i / - {A6) 

J t=constant 

with 5t denoting derivative with respect to time t . 

The in— and out— states of $ are obtained at two extremes of the 
space-time where space-time is asymptotically Minkowskian. These are 
connected through Boglubov transformations 

This scalar product (A6) shows that ($^", = 0, ($^", $|") = 
1 = $1"*). So Prom (A7), it is obtained that 

a, = ($rw,$rw),/?fc = ($rw,$rw) (^8,9) 

obeying the condition 
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i«fei'-i/9r=i. (^10) 

Case of spin-1/2 Dirac fermion 

The equation of Dirac field with mass is given as 

[ij^D^-m^]iP^O {All) 

with I = V^,D^ = + zK^.K)9upKrf Here ^ = KTia = 
0,1,2,3), p''^ = Khlr]"'^ {i^.Y] = 2^'"" and = 27}"^ with 

7'^(7") being Dirac matrices in curved(flat) space-time. 

Further are decomposed in discrete modes k and spin (s/2) with 
s = ±1 as 

oo 

Mlk,s + C?Lfe,>77fe,.] (^12) 

s=±l k=—oo 

oo 

= E E [^ik,sl%ls + i^nk,sfd.u,s] (A13) 

s=±l fc=— oo 

with hk^s{b\^s) and d_k,s{d^-k,s) being annihilation (creation) operators 
for positive (negative) energy particles respectively. Further, using 

V'/M = /fc,.(r/)e~'^-X, (^14) 
= gkAvy''-^: (Alb) 

where 

















/0\ 









1 












Ml = 












1 


,U-l = 






















Here also, in- and oiii- states of ijj are obtained at two extremes of 
the space-time where space-time is asymptotically Minkowskian. 
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The Bogolubov transformations for arc given as 

^-k-s — "fe,sV"jfc,sJ Pk,s^-k-s- 

The Bogolubov coefficients ctfe^^ and j3k,s satisfy the condition 



where 



and 



o^k,s = / dMTk,s^Z'l (^17) 

^ t=constant 



J t=constant 
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